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Abstract 

The conversion of second-class constraints into first-class constraints is used to 
extend the coordinate-free path integral quantization, achieved by a flat-space Brow- 
nian motion regularization of the coherent-state path integral measure, to systems 
with second-class constraints 



(N ■ 1 Second-class constraints 

O, 

On performing the Legendre transformation for generalized velocities of a Lagrangian 

of a dynamical system, one very often gets relations between canonical coordinates and 

momenta that do not involve time derivatives. They are therefore not equations of motion 

and are called primary constraints [|1[] . The primary constraints should be satisfied as time 

proceeds, which leads to further conditions on dynamical variables known as secondary 

constraints [|T|. 

Let ip a = <f a {0) — be all independent constraints (primary and secondary) in the 

system; here 9 % , i = 1,2, ...,2N, denote canonical variables that span a Euclidean phase 

space of the system. The canonical symplectic structure is assumed on the phase space 

{9\9^} =uj % i\ one can, for instance, set q n = 9 2n ~ 1 and p n = 9 2n , n = 1,2, ...,N for 

canonical coordinates and their momenta, then {p n , q m } = 5™ and other components of 

the canonical symplectic structure are zero. Let H(9) be the canonical Hamiltonian of 

the system. Since ip a is a complete set of constraints, 

Va = { V a,H}=C b a (9)Vb~0, (1.1) 

where the symbol ~ implies the weak equality that is valid on the constraint surface 

'-Pa = 0. 

Systems with constraints admit a more general dynamical description where the Hamil- 
tonian can be replaced by a generalized one H? = H+\ a (9, t)ip a (0) with A a being arbitrary 
functions of 9 % and time: 

9 l = {9\ H T } w {9\ H} + X a {9\ Va } . (1.2) 
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The condition ( |1.1| ) with H replaced by Ht enforces some restrictions on A a . Indeed 

0a « {(Pa,^b}^ b = , (1.3) 

that is, the number of independent functions A a is determined by the rank of the matrix 
{(p a , tp b } = A ab on the surface of constraints. If the rank is zero, i.e. {ip a , <p b } = C^ b ip c ~ 
then all the A a are arbitrary, and the solutions 9 l (t) to the equations of motion ( |1.2|) would 
contain arbitrary functions A a . Such constraints are called first-class constraints |l[ and 
transformations of 9 t (t) generated by A a — > A a + 5X a are known as gauge transformations. 
The latter implies in particular that the dynamical system has more non-physical canon- 
ical variables in addition to those fixed by the constraints (p a = 0. They can be removed 
by specifying A a (by gauge fixing). 

If the rank of A ab on the surface (p a = is equal to the number of independent 
constraints then det A ab ^ and all the A a must be zero. The system admits no gauge 
arbitrariness. Such constraints are with an arbitrarily large diffusion constant called 
second-class constraints jjfl. Of course there could well be a "mixed" case when the 
system possesses second and first class constraints (the matrix A ab is degenerate, but 
A a b^0). In what follows only second-class constraints are considered. 

The number of second-class constraints must be even because the determinant of an 
antisymmetric matrix of odd order is always zero. So we set a = 1, 2, ... , 2M, that is the 
system has only 2(N — M) physical canonical variables that describe dynamics on the 
physical phase space determined by (p a (9) = 0. 

2 Local parametrizations of the physical phase space 

One can introduce a local parametrization of the constraint surface 9 l = 9 l { r d) where 
i? a , a = 1, 2, ... , 2(N — M), are chosen so that ip a (9(i))) identically vanish for all values 
of d a . The local coordinates d a span the physical phase space. They may also serve 
as local symplectic variables on the physical phase space, provided there is an induced 
symplectic structure on it. To obtain an induced symplectic structure for an ordinary 
change of variables on phase space, one would have to invert the relations 9 % = 9 % ($) 
and calculate the Poisson bracket of d a . This procedure is incorrect for systems with 
constraints. The inverse relations d a = i9 a (8) are determined modulo the constraints 
ip a = 0. Since {9\<p a } ^ 0, the Poisson bracket {{} a ,'d 13 } is ambiguous and does not 
induce the right symplectic structure. Recall that due to the same reason ({9\(p a } ^ 0), 
the constraints should not be solved before calculating the Poisson bracket in equations 
of motion (|1.2j ). 

The problem is resolved by means of the Dirac bracket that reads as [|T| 

{A, B} D = {A, B} - {A, ip a }A ab { Vb , B} (2.1) 

for any A and B, and A a f,A 6c = 5^. The Dirac bracket possesses three important properties 
|T|. First, it satisfies the Jacobi identity, the Leibnitz rule and is antisymmetric, therefore 
it determines a symplectic structure 

{9\ 9 j } D = ^- {9\ <p a }A ab {<p b , 9 j } = u%{9) . (2.2) 



Second, it vanishes for any A{6) and any of the constraints <p a , 

{A,<p a } D = 0. (2.3) 

As a consequence of ( |2.3| ) we deduce the third property 

A = {A,H} d k{A,H}, (2.4) 

i.e. the equations of motion are not affected by the replacement of the Poisson bracket 
by the Dirac one. 

The advantage of using the Dirac bracket is that one can solve the constraints at any 
stage of calculation, before or after calculating the brackets. The latter is guaranteed 
by (|2.3|) . In particular, given a set of local parameters i9 a = d a {9) spanning the surface 
ip a = 0, the induced symplectic structure is unique 

{& a ,^} D = u afi {-&). (2.5) 



We remark that the Dirac symplectic structure (|2.2j ) is degenerate. Its non-degenerate 
part is given by ( |2.5|) in local coordinates d a . There is an infinite number of choices of local 
symplectic coordinates on the physical phase space. One can find such a parametrization 
for which the induced symplectic structure (|2.5|) has the canonical form u a ^ =uj a/3 . 
The latter follows from the Darboux theorem, and the corresponding $ Q are Darboux 



coordinates for the symplectic structure (2.5). But even the Darboux coordinates are not 



unique, since they are determined up to a canonical transformation. 

Thus, classical dynamics of second-class constrained systems exhibits a physical-phase- 
space reparametrization invariance. 

3 Examples of second-class constraints 

Consider a Lagrangian of the form 

L = lx 2 + yF(x), (3.1) 

where x is a radius vector in 1R .It describes a motion of a point-like particle of unit mass 
on an N — 1- dimensional surface determined by the equation -F(x) = 0, as follows from 
the Euler-Lagrange equation for the Lagrange multiplier y. The canonical momentum for 
y vanishes yielding the primary constraint 

dL 
Doing the Legendre transformation for x, we arrive at the Hamiltonian 



H = -p 2 - yF(x) + X Py , (3.3) 



where p is the canonical momentum for x, and A is an arbitrary function of canonical 
variables and time. Its occurrence in Q3.3|) is due to ( |3.2|) (since p y vanishes identically, 



the canonical Hamiltonian is determined up to any function that vanishes as p y — 0). 

To find secondary constraints, one should calculate the Poisson bracket of the primary 
constraint (13.21) with the Hamiltonian (137 



p y = {p y ,H}^F(x) = 0, (3.4) 

that is, the motion is indeed constrained to the surface F — 0. We should continue 
checking the dynamical self-consistency for the constraint 



F(x) = {F, H} « (p, dF) = . (3.5) 

Equation ( |375| ) determines a new constraint. Let us denote the constraints Q3.2p, ( |3.4[ ) 
and (|3.5|) as (^1,2,3, respectively. Then there must be 

03 = Ws, H} « PiPjdityF + y{dFf = (3.6) 



The new constraint ( |3.6j ) is denoted by ^4. The theory does not have more constraints 
because the condition 04 = {^4, H} « yields an equation for an arbitrary function A, 
rather than for canonical variables. It is not hard to be convinced that all four independent 
constraints ip a are of the second class, i.e. det A ab = det{(p a , (Pb}^ 0. 

A geometrical meaning of (|3.4j) and (|3.5|) is transparent. Equation (|3.4|) implies that 



the particle moves along the surface F — 0. Equation ( |3.5| ) means that the particle 
momentum remains tangent to the surface F = during the motion. Note that the 
vector dF is locally transverse to the surface F = 0. The constraints ( |3.2j ) and ( |3.6| ) are 
artifacts of constructing the Hamiltonian formalism from the Lagrangian (|3.1|) where y is 
not a dynamical variable, rather it is a Lagrange multiplier used to enforce the constraint 
(|3.4j ) in the Lagrangian formalism. Since in the Legendre transformation the variable y 
is treated as an independent dynamical variable, the associated Hamiltonian formalism 
exhibit two extra constraints ( p.2|) and (|3.6j ) to suppress dynamics of y and p y . 

In fact, we may start right from the Hamiltonian formalism ( |3.3| ), setting in it A = 
and y to be an arbitrary function of p and x. Then the constraint ( ft.4| ) should be regarded 
as the "primary" constraint. In this simplified approach equation ( |3.6j ) is not a constraint, 
but the equation for an arbitrary function y (the Lagrange multiplier). The Dirac bracket 
formalism leads to the same answer for the symplectic structure on the physical phase 
space, so we prefer the simplified Hamiltonian formalism. In general, given a set of the 
second class constraints ip a to be imposed on the motion generated by the Hamiltonian H s 
of a system under consideration, one can consider a generalized Hamiltonian formalism 

H = H s (e) + \ a {6,t)ip a (6) . (3.7) 

The consistency conditions (p a = {^> a , H} w yield equations for arbitrary functions A a 
whose solutions are 

\ a = A ab {p b ,H s } . (3.8) 



For instance, to describe the motion along a surface F = in the Hamiltonian formalism, 
one can take H s = p 2 /2 and two constraints ipi = F and ip2 — (p 5 9F) and consider the 
generalized Hamiltonian dynamics (|3.7|) . 

A symplectic structure on the physical phase space determined by the constraints 
(fi 7 2 = is induced by the Dirac bracket 

{xi,Xj} D = ; (3.9) 

{xi,Pj}D = Sij-niTij] (3.10) 

{Pi,Pj}D = p k (njdkni - riidkUj) , (3.11) 

where rij = rij(x) = diF/\dF\ is a unit vector that coincides with the normal to the 
surface when x is on the surface. The symplectic structure is by construction degenerate. 
To construct the induced (non-degenerate) symplectic structure on the physical phase 
space, one should introduce a local parametrization of the constraint surface cpi t 2 = 
and calculate the Dirac bracket for local coordinates spanning the constraint surface. 
For instance, given a local parametrization of the surface F(x) = in the form x = 
x F (w), u G IR , the physical momenta are pi = e^(u)p Ua where the vectors e a (u) form 
a basis in the tangent space of the surface. In particular, one can take ef(u) = dxf /du a . 
From the identity F(x F ) = follows the orthogonality relation (n(x F ),e a ) = 0. The 
variables p Ua ,u a serve as local coordinates on the physical phase space. The corresponding 
symplectic structure is induced by ( |3.9|) -( j37Ll] ). Another possibility would be to solve 



F = with respect to, say, x^ and tp 2 = with respect to any of the momenta, say, p^, 
i.e. u a = x a ,p Ua =p a , a = 1,2,..., JV- 1. 

Let us illustrate the procedure with the two-dimensional case, the motion on a plane 
constrained to a curve. We introduce the following parametrization of the constraint 
surface 

x = f(«) , p =p u dj , (3.12) 

where u is a parameter on the curve, note that dF ~ Td u f along the curve F = 0, where 
Tij = Eij = —Eji, E\2 = 1. The variables u and p u are local coordinates on the physical 
phase space. We get 

Aafc = Wa^ b }=T ab (dF) 2 , (3.13) 

A ab = -T ab {dFY 2 . (3.14) 

Choosing some function u = m(x) (e.g. one can simply invert the relation x\ = fi(u)) we 
obtain 

p n = 7 (x)(p,mF), (3.15) 

where 7(x) depends on the choice of w(x). Hence 

{u,p u } D = 7{x){du{x),TdF)\ F=0 = T(u) . (3.16) 

The Darboux transformation for the physical symplectic form reads {u,p u ) — > (u,p u /T). 
In general, the constraint surface may have a nontrivial topology, which must be taken 
into account when studying the dynamics in local symplectic coordinates. Consider, for 



example, a particle on a circle, -F(x) = x 2 — R 2 . The radial motion is frozen and the other 
constraint is (p, x) = 0. A natural parametrization of the physical phase space (according 
to ( [Tig) ) is 

x = Re ujT X , p = p u Tx(w)/R (3.17) 

where \i — $ii an d w = tan _1 x 2 /xi, p u = (p, Tx). Calculating the Dirac bracket we 
obtain 

{u,p u } D = l, (3.18) 

i.e. uj and p w are Darboux coordinates on the physical phase space which has the topology 
of a cylinder because u is a cyclic (compact) variable, uj G [— tt, tt). One can also choose 
u — X\ and p u = pi as local coordinates on the physical phase space . The corresponding 
symplectic structure assumes the form 

u 2 
{u,p u } D = l-—. (3.19) 



The Darboux transformation for ( |3.19|) is 



u = Rcoslj , p u = ——pujsmu . (3.20) 

R 

Now the topology of the physical phase space is hidden in singularities of the symplectic 
structure: It vanishes at u — ±R. 

The examples illustrates an arbitrariness in choosing a parametrization of the physical 
phase space: The parametrization is determined up to a general coordinate transformation 
on the physical phase space. The Dirac bracket ensures a covariance of the Hamiltonian 
dynamics with respect to such transformations. This covariance of the classical dynamics 
is lost upon canonical quantization as we now proceed to demonstrate. 

4 Ambiguities in the canonical quantization of second-class con- 
straints 

The canonical quantization of a classical system implies that the canonical variables 6 l 
become hermitian operators 9 l that act in a Hilbert space and obey the canonical com- 
mutation relations [8 l ,6^] = ih{8 l ,8i} = ih u % K The recipe is generally correct only in 
Cartesian coordinates 0. Though the canonical variables Qi have been assumed to be 
Cartesian, the quantization postulate should be modified when second-class constraints 
are present. The point is that the conditions tp a = tp a {9) — can not be imposed 
on the operator level because they would be in conflict with the commutation relations 
[ip a , 0b] 7^ 0: The constraints cannot be solved before calculating the commutation rela- 
tions. The problem is resolved by replacing the Poisson bracket by the Dirac one in the 
canonical quantization postulate, that is [||] 

[8\ 9 j ] = i%{6\ 9 j } D \ 9=§ = ihJi0) . (4.1) 



Since the Dirac symplectic structure is degenerate, not every canonical variable can be 
made an operator (e.g. if <pi = q — 0, ip 2 = p = 0, then the Dirac rule leads to [q,p] = 0, 
i.e. the canonical variables are commuatative and, therefore, remain c-numbers upon 
quantization). For a generic second-class constrained system, the Dirac commutation 
relations ( |4.1| ) are constructed so that the operators of constraints commute with canonical 
variables [tp a , 9 l ] = and, hence, can be given any c-number value, in paticular, ip a = 0, 
enforcing the constraints on the quantum level. This comprises the geometrical meaning 
of the Dirac approach. 

The recipe ( |4.1|) is not however free of ambiguities either. The Dirac symplectic 
structure depends on the canonical variables and therefore the replacement 9 l by the 
corresponding operators usually leads to the operator ordering ambiguity. 

An incorrect operator ordering in the right-hand side of (|4.1|) can break the associa- 
tivity of the operator algebra (|4.1|) (the Jacobi identity is violated upon quantization). 
To restore the associativity, terms of higher orders of % should be added to the right-hand 
side of ( f4.1|) . Even after the associativity problem has been resolved in some way, one 
needs still to verify that this solution has not violated the quantization consistency condi- 
tions [9 l ,ip a ] = 0. The latter would involve solving the operator ordering problem in the 
constraints in a way compatible with the operator ordering in the symplectic structure. In 
general this program may be very involved. After all the consistency problems have been 
resolved, one should face a not less difficult problem of constructing a representation of 
the algebra ( [4.1| ) in order to be able to calculate amplitudes (e.g. the evolution operator 



kernel) . 

As an alternative approach one can consider quantization after solving the constraints, 
meaning that the physical symplectic structure ( |2.5|) is to be quantized. It should be noted 
that such an approach, though it resolves the operator ordering in the constraints, still has 
this problem in the symplectic structure ( |2.5| ). Going over to the Darboux variables does 
not help in this regard because canonical quantization is not generally correct in curvilin- 
ear coordinates as has been mentioned before. The problem appears even more serious if 
one notices that the Darboux coordinates are determined up to a canonical transforma- 
tion, whereas canonical quantization and canonical transformations are not commutative 
operations. Thus, such a reduced phase-space quantization is coordinate dependent, and 
in this regard cannot be considered as a self-consistent quantization scheme. 

As an illustration, consider the canonical quantization of the Dirac symplectic struc- 
ture for a particle on a circle. Here n = x/|x|, and canonical quantization of the Dirac 
bracket ( |3.9|) - ( 3.11 ) yields the following commutation relations 



[xuxj] = 0, (4.2) 

[xj,p k ] = ih(5 jk 72~), (4.3) 

X 

\Pj,Pk) = ih—ipjXk - p k Xj) . (4.4) 



The operator ordering problem appears only in (|4.4| ) and can be resolved by putting all 
the pi either to the left or to the right of the X{ in the right-hand side of ( |4.4|) . Note that 



x commutes with all the canonical operators, so it does not matter where it is placed 
in the right-hand side of (|4.4j ). One can check that this ordering is compatible with the 
hermiticity of the canonical operators and provides the Jacobi identity (associativity) in 
quantum theory. In this particular model, the operator ordering problem in the constraints 
is not of any relevance. Indeed, (p\ — x — R 2 does not have any, and commutes with all 
the canonical operators. According to (f4.3| ), the ordering correction to tp 2 is a c-number 
and, hence, does not affect the relations [v?2,%] = [02, Pj] — 0. 

The algebra (|4.2|) - (f4.4|) has a representation in a space of 27r-periodic functions ip(u + 
2tt) = ip(uj) 

Xiip{uj) = Rcosu}ip(ij) , %2ip{uj) = R sin L)ip(uj) ; (4.5) 

ih 

Piip(u) = — (sin udu + du sin u)i()(uj) ; (4.6) 

ZJrC 

ih 
p 2 ip(u) = -— (cos ud^ + d^ cos uj)ijj(uj) ; (4.7) 

(V1IV2) = rdur^2- (4.8) 

Jo 

In this representation (p\ = and 02 is a c-number determined by the chosen operator 
ordering. The physical Hamiltonian assumes the form 

1 h 2 

H ph = 2]pPl + g^2 ' P» = ~ iW " ■ ( 4 - 9 ) 

In addition to the kinetic energy opertor on the circle, the physical Hamiltonian contains 
a "quantum" potential that has occured through the Dirac degenerate commutation re- 
lations. For a generic manifold, the Dirac approach leads to a quantum potential that 
depends on position on the manifold. 

A similar quantum potential was also predicted in the framework of the path integral 
qunatization on manifolds [0], and found to be proportional to the scalar curvature of 
the manifold. It is interesting to observe that canonical quantization of the Dirac bracket 
leads to a different prediction. For an iV-dimensional sphere, the scalar curvature potential 
reads ah N(N — 1)/R 2 with a being a constant, i.e., it vanishes for a circle, whereas the 
embedding of the iV-sphere into 1R +1 and canonical quantization of the Dirac bracket 



( |4.2j ) - ( |4.4| ) would yield the other form of the vacuum energy h N /8R . Note that the 



algebra (|4.2j ) - (|4.4j ) applies to quantum motion on the iV-sphere as has been shown above. 



Its representation is easy to find by going over to the spherical coordinates and thereby 
obtaining the physical Hamiltonian. 

Let us turn to quantization in Darboux variables. The canonical quantization of the 
physical symplectic structure ( |2~5"1) in the Darboux coordinates would lead to a different 
Hamiltonian 

Kh = ^fpPl ■ (4-10) 

In this case the "extra" quantum potential is not unique at all. The source of troubles 
is that the Darboux variables (or any set of canonical coordinates parametrizing the 
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physical phase space) are determined only up to a general canonical transformation and, 
hence, are generally associated with non-Cartesian coordinates for which the canonical 
quantization is not generally consistent @. For instance, one can choose an alternative 
parametrization of the physical phase space by going over to new canonical coordinates 
(u,pj) — > (u — sinu;, p u = p w / costu), {u,p u } = 1. After canonical quantization, [u,p u ] = 
ih, the physical Hamiltonian would exhibit an operator ordering problem which has no 
unique solution. Given a particular operator ordering in the Hamiltonian, the change of 
variables u = sinuj would not generally lead to (|4.10 ), rather the quantum Hamiltonian 



will have an uj dependent potential of order % . This quantum potential is fully determined 
by the physical phase space parametrization chosen. In this regard quantization of the 
Dirac degenerate symplectic structure Q2.2|) in a flat phase space looks more preferable 
because different parametrizations of the physical phase space are associated with different 
realizations of the same algebra of commutation relations for the canonical variables. 

We would also like to mention that we do not consider the physical origin of the second- 
class constraints and assume the latter to be given a priori. In general, the applicability 
of the Dirac formalism to concrete physical sysytems can be questioned. For instance, the 
motion on a sphere can be physically interpreted as a motion in a thin spherical layer so 
that the radial motion is confined by a spherically symmetrical potential well. By going 
over to spherical coordinates in the free-particle Schrodinder equation, one can show that 
in the limit when the layer thickness is much less than the sphere radius, the quantum 
potential is given by the quantum centrifugal barrier h N(N — 2)/8R 2 acting on a particle 
in the thin spherical layer, which is different from the one predicted by the Dirac or path 
integral approaches. 

So, it is an open question whether or not the Dirac quantization scheme for second- 
class constraints can be applied to a particular dynamical system. 

5 Quantization of second class constraints via the abelian con- 
version method 

We have seen that the quantization of the Dirac bracket poses a few problems: the operator 
ordering in the commutation relation algebra (the associativity problem), the problem of 
finding a representation of the Dirac commutation relations and the ordering problem in 
the operators of constraints. Quantum dynamics depends on a particular solution to them 
and, generally, is not unique. The first two problems are most difficult. So one should 
develop a formalism that would allow one to avoid them. Such a formalism is known 
as the conversion of second-class constraints into first-class constraints by extending the 
original phase space by extra (gauge) degrees of freedom Q. 

Let us extend the original phase space of a system with 2M independent second-class 
constraints by adding to it 2M independent variables (f) a with the canonical symplectic 
structure 

{0 a ,0 b }=c> b , {0 a ,^} = O. (5.1) 

The original second-class constraints ip a {@) are then converted into abelian first-class con- 



straints a a = a a (8, <fi). Their explicit form is determined by a system of first-order differ- 
ential equations 

K,C7 6 } = (5.2) 

with the initial condition 

<Ta(e,<f>=0)=<Pa(9). (5.3) 

A dynamical equivalence of the original second-class constrained system to the abelian 
gauge system is achieved by a specific extension of the original Hamiltonian of the system, 
H s (6) -> H 8 (9,<j)), such that 

{H s ,a a } = 0, H S (6,<P = O)=H S (0) . (5.4) 

One can show that equations of motion generated by the extended action 

S = Jdt (hfi u ah ^ h + ~0* UijP -H s - AV a ) , (5.5) 

are equivalent to those generated by the original action 

s = Jdt(~e i u ij e j --H s -\ a pa) ■ (5.6) 

We remark that in contrast to the theory ( |5.6| ), the Lagrange multipliers A a in the gauge 
theory ( |5.5| ) are not determined by the equations of motion ( the matrix A ab = {a a , a b } 
vanishes according to ( |5.2j )). A choice of A a implies gauge fixing. In particular, one can 
always choose A a so that <p a = on the constraint surface a a = for all moments of time. 



With this choice the equations of motion of the system (5JJ) become the equations of 
motion of the original system. 



Equations (|5.2| ) and (|5.4j ) are not easy to solve for generic (f a and H s . However, if at 



least one set of Darboux variables for the Dirac bracket is known, then the solution can 
be found explicitly [f| . In particular, for a particle on a circle we find 

a x = <p 1 + P = x 2 -R 2 + P , (5.7) 

a 2 = <p 2 + 2x 2 Q = (x, p) + 2x 2 Q , (5.8) 



where 1 = Q, <ft 2 = P and {Q, P} = 1. To find a solution to ( |5.4|) , one can make use of a 
simple observation that 0\p, and p^ = (p, Tx) could be regarded as canonical momenta for 
Q, hir/R and uj = tan -1 x 2 /xi, respectively. In these new canonical variables equation 
J) is greatly simplified and we obtain 

1/ a\ (p,Txn 

Hs ~ 2 {^TW + -^TW) ■ (5 ' 9) 

When P = Q = 0, o\ + R 2 = x 2 and a 2 = (p, x) so that H s turns into the free-particle 
Hamiltonian H s = p 2 /2 (recall that the vectors x and Tx form an orthogonal basis on a 
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plane). We remark also that formulas (|5.7|) - (|5.9| ) apply to an iV-dimensional rotator (a 
motion on the iV-dimensional sphere); one should only replace the squared total angular 
momentum of the rotator by its N- dimensional analog: (p, Tx) 2 — > L 2 = J2a{P^T a x.) 2 
where T a are N x N real antisymmetric matrices, generators of SO(N). 

The first-class constraints ( |5.7| ) and ( p.8| ) generate gauge transformations on the ex- 
tended phase space which are translations of Q and dilatation of the radial variables r, 
while the angular variables and its momenta remain invariant. Thus, Q and r are pure 
gauge degrees of freedom, and the angular variables comprise guage invariant physical 
degrees of freedom as expected. 

The gauge transformations are also canonical transformations with generators being 



01,2- An operator of finite canonical transformations generated by the constraints ( |5.7|) 
and ( |5.8| ) has the form 

exp(£ a ada a ) , ada a = {cr a , •} . (5.10) 

Applying it to the canonical variables in the extended phase space, one finds 

Q -» Q~Zi = Qt, P^P+(l-e- 2 « 2 )x 2 = P € ; (5.11) 

x -► e^ 2 x = x c , p^e 6 p + 2x(Qe 6 -(Q-6)e _6 ) = p € • (5.12) 



From ( p.ll| ) follows that there always exists a choice of the gauge parameters £ a such that 



Q = P = for all moments of time. 

After the conversion has been made, the system can be quantized according to the 
Dirac method for the gauge theories. Namely, all the canonical variables of the extended 
Euclidean phase space become operators obeying the standard Heisenberg commutation 
relations 

[fa, § k ] = ih ^ jk , [<T, (j> b ] = ihu ab , (5.13) 

while the operators of constraints a a select physical states 

cr a y ph = 0. (5.14) 



Equation ( |5.14j ) means that the physical states must be invariant under gauge transfor- 



mations generated by first-class constraints. In particular, for the rotator we find that so- 
lutions to the Dirac constraint equations fl5.14|) are given by functions ^ p h = f(Q, r)if)(u)) 
where f(Q, r) is uniquely fixed by Q5.14J ), while ip(u>) is an arbitrary function of the polar 
angle on a plane. For the iV-dimensional sphere, the physical Hilbert space consists of 
functions on the sphere. So, it is the Hilbert space of a quantum rotator as expected. 

Thus, the problem of quantization of second class constraints has a natural gemetrical 
solution in the framework of the conversion method. The technical difficulties do not 
disappear completely; they are now associated with solving the conversion equations (|5.2|) 
and ( |5.4j ). Nevertheless, the approach may be simpler than the original Dirac approach. 
Even in the case of the rotator, the representation problem for the algebra ( f4.2|) -( f4~4|) 
does not appear to be a feasible task if the geometrical origin of this algebra is unknown. 

An important advantage of the conversion method is that it does not rely on any 
particular parametrization of the physical phase space. For this reason we shall adopt it 
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as starting point to develop a path integral formalism for second-class constrained systems 
which is invariant with respect to the parametrization choice of the physical phase space 
and, in this sense, to achieve coordinate independence of the quantum theory. 

6 The projection method 

We assume the operators a a to be hermitian and that they generate unitary transforma- 
tions in the total Hilbert space. Since by construction they commute with the Hamil- 
tonian, the total Hilbert space of an abelian gauge system obtained by the conversion 
procedure can always be split into an orthogonal sum of a subspace formed by gauge 
invariant states ( pM4| ) and a subspace that consists of gauge variant states. Therefore an 
averaging of any state, being a linear combination of eigenstates of the Hamiltonian, over 
the abelian gauge group automatically leads to a projection operator onto the physical 
subspace of gauge invariant states: 

in a 



V = / 5Me' lU Ua (6.1) 

where 5 a Q is a normalized measure on the space of gauge transformation parameters. If 
the spectrum of the constraint operators a a is not discrete, then the parameters Q a range 
over a non-compact domain. In this case we adopt a certain regularization of the measure 
5 a Q which provides (see also Section 9) 



/ 



5 a Q = 1 (6.2) 



and, hence, V is the projection operator V 2 = V such that 

V% h = % h , Vy nph = (6.3) 

for any gauge invariant state $ p h and any gauge variant state ^ np h (by definition, a a ^nph 7^ 
0). Its kernel is determined as the gauge group average of the unit operator kernel 

9',(j)') ph = {e\(f)"\V\e\<p') = f 8 a n{6", <p"\e- inafTa \9',<t>') , (6.4) 



where \6,<f)) is the coherent state defined as 

\e, <p) = ex P (ie j u jk e k + %<\> a u ab ft) |o> (6.5) 

with |0) being the ground state of the harmonic oscillator. 

Accordingly, the physical transition amplitude in the coherent-state representation is 
obtained from the unconstrained one by averaging the latter over the gauge group 

(9",(f)",T\6',(l)') ph = f 5M(0", <P" ,T\e- iQa&a |0',0') (6.6) 

j^ m (e",^,T\9, ( p)(9, ( p\V\9', ( p') . (6.7) 
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The unconstrained transition amplitude is given by the coherent-state path integral 

T 



T\e',(f)') = fv6V(f)exp lifdt 

V o 



911 Jjl rr\nl J.' 
,0 , rT "' 



h% + ^<p a <Pa-h s (e,. 



(6.8) 



with the boundary conditions 9(0) = 9', 9(T) = 9" and (f)(0) = <ft', 4>(T) = <ft"; here 
9i =u> i j9 : * , (p a =^>ab(p b and h s is the lower symbol for the operator H s 

^ = / (2^(2^ ^' 0)l ^ 0)(0 ^ 1 - (6 ' 9) 

Dividing the time interval T into n pieces e = T/n and taking a convolution of n kernels 
( |6.7j ) where T — > e, we arrive at the coherent-state path integral representation of the 
physical transition amplitude 

(9", <P", T\9', <\jf h = f V9V(t)VC{uj)e l ^ dt ^ e%+ ^ a ^- uaaa -~ ha ) . (6.10) 

The measure T>C(uj) for gauge variables, being the product of the local measures 5 a uj a (t), 
provides the projection at each moment of time. The action in the exponential in (|6.10|) 
coincides with the classical action ( |5.6|) up to possible operator ordering terms H s — h s = 
0(h). It is invariant with respect to gauge transformations generated by a a 

59 l = C^a a 9 { , 5<f) a = i\da b (f) a , 5u a = -£ a , (6.11) 

where the infinitesimal functions of time £ a satisfy zero boundary conditions 

£ a (0) = C(T) = , (6.12) 

which ensure that the boundary terms occurring upon varying S vanish. 

To obtain the corresponding path integral on the physical phase space, one usually 
has to integrate out all the gauge variables u and non-physical degrees of freedom. For- 
mally, it can be achieved by going over to new canonical variables such that the abelian 
constraints a a become new canonical momenta. Due to the gauge invariance the Hamil- 
tonian h s is independent of the corresponding canonical coordinates. Since the Liouville 
measure constituting the formal path integral measure is invariant under canonical trans- 
formations, the integration over non-physical variables becomes trivial. Yet, the gauge 
transformations in the new variables are translations of canonical coordinates for the new 
momenta a a , so the gauge average can also be done explicitly. Having done this, one 
seems to obtain a path integral over the physical phase space parametrized by a certain 
set of canonical variables. Note that the canonical transformation discussed above is not 
unique and is determined only up to a general canonical transformation on the physical 
phase space. On the other hand, we have seen in section 4 that quantum theory may well 
depend on a particular parametrization of the physical phase space, which is in conflict 
with the formal coordinate invariance of the path integral measure. 
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To explain the contradiction we observe that the above procedure of reducing the 
path integral measure onto the physical phase space relies on the formal invariance of the 
conventional Liouville measure with respect to canonical transformation. Unfortunately, 
this is a wrong assumption. A typical example is a Hamiltonian dynamics on a phase-space 
plane. By a canonical transformation one can always locally turn the Hamiltonian into a 
free particle Hamiltonian. So assuming a formal invariance of the path integral measure 
with respect to general canonical transformations we would arrive to a contradiction 
that every quantum dynamical system with one degree of freedom is equivalent to a free 
particle. 

Thus, in order to obtain a path integral on the physical phase space, the path integral 
measure in the amplitude ( |6.10| ) should be regularized in a way that provides a true 
covariance of the path integral ( |6.10| ) with respect to general canonical transformation. 

7 The Wiener measure regularized path integral 

For Hamiltonian systems without constraints a regularization of the path integral measure 
can be achieved by replacing the conventional Liouville measure by a pinned Wiener 
measure on continuous phase-space paths. The Wiener measure regularized phase space 
path integral for a general phase function G(p, q) is then given by [[JJ 

lim M v \ exp{iJ^[pjq J + G(p, q) - h(p, q)] dt} 

X exp{-(l/2v)jT[p 2 + q 2 } dt} VpVq 

= \\m(2n) N e N » T / 2 J exp{tQp 3 dq^ + dG(p, q) - h(p, q)dt]} d^ w {p, q) 

= (p",q"\e- WT \p',q'} , (7.1) 

where the last relation involves a coherent state matrix element. Here we use the con- 
vention adopted in Section 1 that g J = O 2 ^ 1 and pj = 9 2 ^ (cf. ( p.5| )). In this expression 
we note that Jpj dq 1 is a stochastic integral, and as such we need to give it a definition. 
As it stands both the Ito (nonanticipating) rule and the Stratonovich (midpoint) rule of 
definition for stochastic integrals yield the same result (since dpj(t)dq k (t) = is a valid 
Ito rule in these coordinates). Under any change of canonical coordinates, we consistently 
will interpret this stochastic integral in the Stratonovich sense because it will then obey 
the ordinary rules of calculus. We also emphasize the covariance of this expression under 
canonical coordinate transformations. In particular, if pdq = pdq + dF(q, q) characterizes 
a canonical transformation from the variables p, q to p, q, then with the Stratonovich rule 
the path integral becomes 

(p",q"\e-* nT \p',q') 
= \im(27T) N e N ^ 2 J exp{i/Jfed? + dG{p, q) - h(p, q)dt}} d^ w (p, q) 

= lim M u / exp{ifJV.<f + G~{p,q) - h{p,q)dt}} 

xexp{-(l/2u)Jo[do-(p,q) 2 /dt 2 ]dt}VpVq, (7.2) 
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where G incorporates both F and G. In this expression we have set da(p, q) 2 = dp 2 + dq 2 , 
namely, the new form of the flat metric in curvilinear phase space coordinates. We 
emphasize that this path integral regularization involves Brownian motion on a flat space 
whatever choice of coordinates is made. Our transformation has also made use of the 
formal - and in this case valid - invariance of the Liouville measure. 

In order to fulfill our program of a coordinate-free path integral representation of 
second-class constrained systems, we have to extend the Wiener measure regularization 
to such systems. 

8 The Wiener measure for second-class constraints 

The regularized measure in the path integral ( |6.10| ) is obtained by the replacement 

V6V(f)VC(uj) -► VC(uj)dfi 9 w (6, 0, u) , (8.1) 

where the gauged Wiener measure d[i 9 w is to be found. The Wiener measure regularization 
of the path integral should not violate gauge invariance, therefore, we impose the condition 

5d/4(0,0,u) = O, (8.2) 



where the operator 5 is determined in ( |6.11| ). Since the Wiener measure provides for 



covariance of the path integral (|6.10| ) relative to canonical transformations, we perform a 
canonical transformation in ( |6.10| ) such that a a become new canonical momenta 

(e\r)^(na = va,y a ,# a ); (8.3) 

here y a are canonical coordinates for ir a and d a are canonical symplectic variables on the 
physical phase space (Darboux variables for the physical symplectic structure (|2.5D ). In 
the new variables the action assumes the form 

S = f (n a dy a + \d a d$ a - u a n a dt - dG - h s di) (8.4) 

and for the Wiener measure we get 

dfi 9 w (6, (j), w) = dfi g w (iT, y, ■&, u) . (8.5) 



The gauge transformations (|6.11 ) leave all the new variables untouched except the y a 
which are shifted 

6ir a = 0, ^ a = 0, 5y a = -C. (8.6) 

By the change of integration variables 

y a {t) - y\t) - f T dt'u a {t') , (8.7) 



one can remove the dependence on u a of the integral ( |6.10|) for all intermediate moments 



of time < t < T. However the initial values of y a are not integration variables and 
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therefore the average over u a does not disappear without a trace. Let us make a change 
of gauge variables u a — > uj a , where the new variables uo a satisfy the boundary condition 
uj a {T) = 0. Note that we are free to add any constant to uj a because it does not affect 
the derivative uj a ] the boundary condition fixes this arbitrariness, providing a one-to-one 
correspondence between the old and new gauge integration variables. With this choice 
equation (|8.7|) assumes a simple form y a (t) — > y a (t) + oj a (t). At the boundary t = we 
have 

y a (0) -> y a (0) + n a , uj a (0) = tt a . (8.8) 

Therefore the average measure for gauge variables is reduced to a single average over fl a , 

VC(u) -»■ 5M (8.9) 

because JVC(u) = 1. 

Equation (|8.2| ) is easy to solve in the new canonical variables 

dfiw(ir,y,$,u) = dj2 9 w (ir,y - lu,$) , (8.10) 



where the gauge variables u have been replaced by their time derivatives as in ( |8.8| ). Note 

that under gauge transformation S(y a — u a ) = £ a — 5u a = because 5u a = —^ a in 

accordance with ( |6.11| ) and the replacement u a —* uj a . Clearly, the further transformation 



of the y-integral to the new variables y — u removes the dependence of the Wiener measure 
on the gauge variables for all intermediate moments of time, i.e., 

dp? w {n : y-u,$)^djl 9 w {'K ) y,'&) (8.11) 

in the path integral ( |6.1U| ). As a result of these two canonical transformations the entire 
dependence of the path integral measure on gauge variables is reduced to a single average 
over a gauge orbit of the initial phase-space point with some phase factor determined by 
the phase function / dG of the canonical transformation ( |8.3| ). That is, we have recovered 
the projection formula (|6.7|) where the projection operator kernel is given by 



(9,4>\9',4>r h = SM(9A\0'n,4>n)e lOie '^' n) , (8.12) 



with 0q and 0^ being gauge transformations of the extended phase space variables gen- 
erated by (|5.10| ) with £ = Q and G is G(t = 0) written in the initial canonical variables. 
The Wiener measure regularized path integral ( p.8| ) involved in the projection formula 
( |6.7p should have the flat Wiener measure on the extended phase space according to our 
consideration in the previous section, that is, 

djl 9 w (7i, y, i?) = d/iw(0, 0) . (8.13) 

Having established the relation between djl 9 w and the flat-space Wiener measure d/jJ^, we 
can perform a canonical transformation inverse to (|8.3j) to restore the dependence of the 
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Wiener measure on the gauge variables and thereby to find an explicit form of the desired 
gauged Wiener measure (|8.1| ). Combining (|8.5| ), ( p.lOj ) and (|8.13|) we conclude that 

d\x 9 w {Q, 4>, uo) = d[iw(0u, (f>u>) , (8.14) 

where 

&l = e^^Qi , <pi = e u>baA ' Jb <p a (8.15) 



with uo a = uo a (t). The gauge invariance of the gauged Wiener measure Q8.14D follows from 
the simple observation that 

K = 6<fc = (8.16) 

under the gauge transformation Q6.ll ) where Su> = £ according to the change of gauge 

variables uo — > Co. 

Thus, the Wiener measure regularized path integral for second-class constrained the- 
ories has the form 



r,<p", 



T\e\<p') ph = [vc^fd^ie^^yfoH^^rta-^-K). (817) 

dlfylfimfa) = e( N+M > T / 2 V9V^exp(--±-J Q T ' dt0l + fa)) , (8.18) 



where the limit v — > oo must be taken after calculating the path integral ( |3.17|) . In general 



the gauged Wiener measure (|8.18|) depends not only on Co but also on uo themselves, 



therefore, it is not possible to remove the dependence of the action in (|8.17|) on the time 
derivatives of the gauge variables by changing the gauge variables back Co — > uo in the 
average measure T>C(uo), while maintaining the locality of the gauged Wiener measure 

(ra- 



As an example consider the gauged Wiener measure for the two-dimensional rotator 
(a generalization to the N- dimensional case is trivial as remarked after Eq. ( |5.9| )). The 
canonical transformation (18.3) can be chosen as 



TTi = <j x , y' = Q; (8.19) 

vr 2 = a 2 , y 2 = ln(|x|/it:) ; (8.20) 

d 2 = (p,Tx) , d 1 = tan-^xa/xi) , {^,^ 2 } = 1 . (8.21) 

As expected from a,da a H s = 0, the canonical coordinates y a are cyclic. The gauge trans- 
formations 

y'-ty'-e, y 2 ^y 2 -e (8.22) 

induce gauge transformations of the initial canonical variables (|5.11|) . Setting in (|5.11|) 
C, = uo we obtain the gauged flat metric on the extended phase space that determines the 
Wiener measure 

£ dt(6l + <j>l) = £dt(i>l + ±l + P 2 + Ql) (8.23) 

= [ T dtg AB (A)A A A B , (8.24) 

Jo 
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where A A denotes the set of all canonical and gauge variables (9,cj),u) = (p,x,P,Q,u); 
for the iV-dimensional rotator, p and x are iV-dimensional vectors in (|5.11|) and (|8.24|) . 
Note that the metric gAB depends generally on all the A A , as well as the components gAu 
and g WUJ do not vanish. Thus, the Wiener measure depends on gauge variables and their 
time derivatives. 

Expression (|8.24fj holds for general second-class constrained systems. Its geometrical 



meaning is transparent. The metric qab is, by construction, degenerate along the direc- 
tions traversed by gauge transformations of the A. Hence the gauged Wiener measure 
describes a Brownian motion (with diffusion constant that tends to infinity) in the di- 
rections transverse to the gauge orbits, while the average over the gauge variables with 
the measure T>C(u>) regularizes the path integral along the gauge orbits. An explicit 
construction of this measure is discussed in section 9. 



An unusual feature of the integral ( 3.17 ) is the appearance of the time derivatives of the 
gauge variables in the classical action. This was the price we paid for locality of the Wiener 
measure. One should realize that this is not always the case for the Wiener measure in 
gauge theories. If the canonical transformations generated by first-class constraints were 
linear and preserving a bilinear positive form on the extended phase space, then the 



associated Wiener measure (|8.18|) would have had no dependence on uj but on to only. 



The latter occurs for Yang-Mills type gauge theories ||. In this case the dependence on 
Cj can be removed by a simple change of variables u — > u in the gauge average integral 
without violating the locality of the Wiener measure. 

At first sight, the presence of the time derivatives of the gauge variables in the classical 
action seems to allow for non-physical motion with a = const 7^ (a variation of the action 
relative to u leads to the equation of motion & = rather than just a = 0). One has 
however to bear in mind that equation ( 8.17| ) describes a quantum motion whose gauge 



invariance is ensured by an appropriate average over the gauge variables. As long as the 
measure T>C{iS) provides at least one gauge group average in the time interval < t < T, 
contributions of states with a = const ^ are projected out from the transition amplitude 
in full accordance with the projection formula 



9 The average measure for gauge variables 

The spectrum of the first-class constraint operators that usually occur upon the abelian 
conversion of second-class constraint operators is continuous. Therefore the average mea- 
sure 5 a Q in the projection operator (|6.1|) should be regularized to provide the normaliza- 
tion condition ( |6.2|) . Since the converted constraints are abelian, the projection operator 



( |6.1| ) is the product of the projection operators for each independent abelian generator 
a a . The latter allows us to treat the measure 8 a Q as the product of normalized measures 
for each independent gauge variables Q a . So we can drop the index a and consider the 
measure only for one generator a. 

The gauge transformations are translations of the gauge parameter Q. Hence any reg- 
ularization (a cut-off) of the translation invariant measure dQ would break the translation 
invariance and therefore an explicit gauge invariance of the path integral ( |S.17 ) . In this 
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sense, the regularization would lead to a "gauge-fixing" term in the effective action in the 
integrand in (|8.17|) . The gauge invariance of the amplitude ( fB.17|) is guaranteed as long 



as the regularized measure for gauge variables provides at least one projection onto the 
physical subspace in the time interval t G [0, T]. 

Consider the regularized measure of the following form 



Tfl _rafj2 



S a Q = , — e -~ u dQ, m -> . (9.1) 

V 2n 

Here m is the regularization parameter. Clearly, the measure (|9.1|) is normalized to unity. 
Let \a) be an eigenvector of the generator a. Applying the projector ( p.l\ j to it we find 



V\a) = f 5 a Qe~ ina \a) 

oo 

J dQ e-^ n2 ~^\a) 



oo 

m 




2tt 

— oo 

2 



= e~^\a) . (9.2) 

Taking the limit m — ► in (|9.2j ) we see that for a hermitian operator a, the operator V 
annihilates all eigenvectors of the gauge generator a unless a = 0. In the latter case V 
acts as the unit operator, that is, it is the projector on the physical subspace. 

Adopting the above regularization of the gauge average measure, we replace T in 
the projection formula ( |6.6|) by an infinitesimal time interval e = T/n and construct a 
convolution on n infinitesimal propagators ( |6.6| ). The result has the form ( |6.1(J| ) where 
the gauge variable measure is 



n-l 



VC{uo) = J] J^e-^du, = Afe-fo dt ^ 2 ^ \[ ^ dw{t) , (9.3) 



where u;(0) = Vt (to match the notations in ( |6.6| )). To take the continuum limit we have 



rescaled the gauge variables Uj — > y/eujj with e being the time slicing so that ujj = u>(tj) 
and uj{tj + i) = u>(tj + e). 

To make the gauged Wiener measure a local functional of gauge variables, it was 
proposed in section 8 to change the integration variables u(t) — > u(t). In the time-slice 
approximation of the path integral the change of gauge variables assumes the form 



u 



<j -»• (u j+ i - Uj)/e , j = 0, l,...,n-l . (9.4) 



The "extra" new variable u n = u(T) is fixed by the boundary condition u(T) = as 
suggested in Section 8. In the new variables the measure for gauge variables turns into a 
flat-space Wiener measure for continuous paths pinned at one point 



A/"exp[-^/ dt u 2 j Y[y/m du(t) 



VC m {dj) = AAexp -— / dt Co 2 TT v^ duo(t) , u(T) = , (9.5) 
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where the index m stands to emphasize the dependence of the measure on the regular- 
ization parameter m. With this choice of the measure for gauge variables, we arrive 
at our coordinate-free and mathematically well-defined formulation for the path integral 
representation of the second class constrained systems. 

To conclude the discussion, we note that by construction the limits m —> and v — > oo 
commute in the integral ( 8.17 ). The latter follows from the projection formula ( |B.12| ) to 



which the path integral can be transformed by a change of variables as has been shown in 
Section 8. The amlitude (9, <f)\0' n , <fr' n ) can be calculated at a finite v. For any fixed Q it is 
a regular function of v in the vicinity of v = oo. So, taking its gauge invariant part either 
before the limit v — > oo or after it would yield the same result. It is convenient then to 
make a particular choice of the parameter m to simplify the path integral measure form. 
Namely, we set 

m = l/v , (9.6) 

so that the path integral measure would depend only on one parameter to be taken to 
infinity after performing the sum over paths. Thus, the gauged Wiener measure assumes 
the following (unified) form 

dptty(9 u ,<j> u )VC 1/v (u) = e( N+M > T / 2 e-^fo dt (%+ti+» 2 )veV(j)Vu . (9.7) 

As has been shown in Section 8, by a suitable change of variables one can always remove 



the dependence of the integrand in ( 8.17 ) on the gauge variables u a for all moments of 
time except t — 0. Then the integral over u a yields the kernel of the imaginary time 
transition amplitude for a free motion in the M-dimensional Euclidean space of gauge 
parameters where u a (T) = and o; a (0) = Q a . The integration over the initial values fl a 
weighted with this kernel is precisely the gauge average (|6.6|) regularized as prescribed by 



(9.2). The measure (9/7) describes a Brownian motion with an arbitrarily large diffusion 



constant on the unified space (9, 4>, u), and, in this sense, all degrees of freedom, physical 
and gauge ones, are treated on equal footing in the Wiener measure (|9.7j ). 

We note that Ashworth has also studied Wiener measure regularizations for systems 
with first class constraints @. 
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